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1. INTRODUCTION 

Smart grids are a set of technologies, concepts and approaches, allowing the integration the 
generation, transmission, distribution and use into one internet by full use of advanced sensor measurement 
technology, communications technology, information technology, computer technology, control technology, 
new energy technologies [1]. However, Smart Grid uses digital technology to control grid and choosing the 
best mode of power distribution to reduce energy consumption, reduce costs, increase reliability and also 
increase transparency in the network. Therefore, the system intelligent will have will have a significant 
impact in the fields of finance and economics of the power industry [2]. Although, The traditional network is 
a one-way network in which the electrical energy produced in power plants is channeled to consumers 
without information to create an automated and distributed network of advanced power supplies. 

ED is also applied in the integrated system for scheduling power plants. A few methods have been published 
to solve the ED problem and Optimal Power Flow (OPF). Researchers have published a few methods to solve 
ED and OPF problems. Direct method is accurate and very simple but limited by the quadratic objective 
function [3]. 

The economic dispatch (ED) is one of the power management tools that is used to determine real 
power output of thermal generating units to meet required load demand. The ED results in minimum fuel 
generation cost, minimum transmission power loss while satisfying all units, as well as system 
constraints [4-5]. 

The rise of energy demand and insufficient of energy resources are required for quality and secured 
dispatch [6]. A well-coordinated and optimized power system operation help in satisfying Economic 
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Dispatch (ED) among users of power networks. Hence, studies need to be conducted in order to analyze and 
develop new tools so that the optimization issues in ED could be overcome. Basically, the principal objective 
of load dispatch is to minimize the total fuel cost while satisfying the requirements of some important 
operational parameters. In today’s environment, efficient load dispatch requires not only to schedule the 
power generation at the least cost but also to consider other performance factors to be optimized in power 
flow over the networks. The obligation of social attentions has influenced the reduction of energy 
conservation and pollution emission produced by power plants [7]. 

Facing the electric bulimia experienced by the world and as an urgent and efficient solution is 
sought, it is essential to optimize the cost of producing electricity. As such, tiny cuts costs conceal huge 
potential savings, this is part of this paper, and we look at the overall optimization purposes known economic 
load distribution (OPF) or economic dispatch (ED) [8]. 

The ED is a static problem is to say we must define at a given powers generated by each power plant 
to power a load as economically as possible. To solve this problem the optimization methods are used. 
Conventional optimization techniques [9-10]. have long been applied to solve the ED problem such as 
Quadratic Programming [11-12]. linear programming [13] sequential approach with a matrix framework 
(SAMF) ,[14]. modified Lambda-iteration method [15], Newton Raphson and Lagrangian multiplier (LM) 
algorithms [16], In the real-design cases, the number of decision variables (i.e. power units) of the ED area 
are very large. The objective criterion to be minimized could also have too many local minimum which might 
not lead to the minimum cost and the best generation schedule of power system units. Therefore, efficient 
search algorithms are needed. 

Nature-inspired metaheuristic search algorithms gain a popularity due to their promising 
performance on solving many real-world optimization problems which are complex, nonlinear and multi- 
model. In the past two decades, the literature of metaheuristic search has expanded extensively. 

Some of the well-known metaheuristic approaches are Genetic Algorithms [17], Genetic Programming [18- 
20], Particle Swarm Optimization [21-22], Simulated Annealing [23], Artificial Bee Colony (ABC) [24], 
Cuckoo Search [25-26], etc. 

The remaining organization of this paper is as follows. Section 2 presents the mathematical 
formulation of the ED. Handling of constraints and implementation of the proposed CSA to ED problem are 
addressed in Sect. 3. Section 4 reports results of the proposed CSA method. A number of case studies using 
standard test systems are used to test the proposed method. The comparisons of results between the proposed 
method and existing methods are also carried out in this section. The discussion is followed in Sect. 5. After 
all, the conclusion is given. 


2. MATHEMATICAL FORMULATION OF ED 

The overall research methodology involved in the Economic Dispatch (ED) was classified into four 
stages. The first task was to achieve the objective of the study which was to establish a new technique 
particularly to solve ED optimization problem. The research approach was to design a new optimization 
technique taking some inspiration from the Meta-EP mutation strategy. The development also included 
identifying suitable objective functions which were significant to ED problem along with some constraints as 
discussed in the following section [27]. In order to achieve the research objective, the development of the 
new single objective technique was to be accomplished. The performance of the developed technique was 
evaluated and compared with other techniques namely the AIS and Meta-EP along with Base technique. The 
developed techniques were tested on the standard IEEE 26 and 57 bus system in order to minimize the total 
fuel cost, emission dispersed and system losses. 

The objective function of the ED problem is to minimize the total production cost, which be written 
as: 


MinimizeFy = Y§_, F,(P;) i = 1,2,....,N (1) 
Mathematically, the fuel cost of a thermal generation unit is represented as quadratic function [4]: 

F;(P;) = dj + b;P; + c,P? (2) 

The solution of ED can be highly improved by introducing higher order generator cost functions. 
Cubic cost function displays the actual response of thermal generatorsmore accurately. 

The cubic fuel cost function of a thermal generating unit is represented as follows [28]: 


F,(P;) = a; + bjP, + ¢,P? + d,P? (3) 
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Subject to Real power balance equation The total active power output of generating units must be 
equal to total power load demand plus power loss: 


Lisi Pi = Pp + Py (4) 
where the power loss PL is calculated by the below formulation [4]: 
Py = hy Dk P; BiyP; + Liki Boi P; + Boo (5) 


Generator capacity limits The active power output of generating units must be within the allowed 
limits: 


Pimin < Pi s Pimax (6) 


3. CROW SEARCH ALGORITHM 

In this section, it is explained the results of research and at the same time is given the 
comprehensive discussion. Results can be presented in figures, graphs, tables and others that make the reader 
understand easily [2], [5]. The discussion can be made in several sub-chapters. 

The crow search algorithm (CSA) is a new population-based stochastic search algorithm recently proposed 
by [29]. The CSA is a newly developed optimization technique to solve complex engineering optimization 
problems [30-31]. It is inspired by the intelligent behavior of crows.The principles ofCSA are listed as 
follows [29]: 

a. Crows live in the form of the flock. 

b. Crows memorize the position of their hiding places. 

c. Crows follow each other to commit thievery. 

Crows protect their caches from being pilfered through probability. 

Following the above assumptions, the core mechanism of the CSA consists of three basic phases, 
namely initialization, generate a new position, and updating the memory of crows. At first, the initial 
population of crows represented by n dimension is randomly generated. At iteration t, the position of crow is 
specified by x't = Estee aes »xy and it is assumed that this crow has memorized its best experience thus 
far in its memory mit = [meme jemeaty my; |To generate a new position, crow i select randomly a crow j, 
for example, from the population and attempts to follow it to find the position of its hiding place (mj ). In this 
case, according to a parameter named awareness probability (AP), two states may happen: 

a. State 1: Crow j does not know that crow i is following it. As a result, the crow i will determine the 

hiding place of crow j. 

b. State 2: Crow j knows that crow j is following it. As a result, to protect its cache from being 
pilfered, the crow j will fool crow i by going to another position whitin the search space. 

According to States | and 2, the position of the crows is updated as follows: 


iiter+1 ; iiter iiter _  _i,iter 

xbiter+1 _ ae +1; X fi x (m xe), (7) 
rj = AP!"**'A random positonofsearch spaceotherwise 
Where rj is a uniformly distributed fuzzy number from [0; 1] and AP! teT denotes the awareness probability 
of crow j at iteration iter. Finally, the crows update their memory as follows: 

miitertt — io if fxr )is better than f(m! er ) (8) 

m?ter otherwise 

Where f(-) denotes the objective function value. It is seen that if the fitness function value of the new position 
of a crow is better than the fitness function value of the memorized position, the crow updates its memory by 
the new position. The above process is repeated until a given termination criterion (itermax ) is met. Finally, 
the best solution of the memories is returned as the optimal solution found by the CSA. The main steps of the 
CSA are outlined in Algorithm: 
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Algorithm Crow Search Algorithm 
i: Randomly initialize the position of a flock of (Np ) crows in the search 
space. 

2 Evaluate the position of the Crows 
3 Initialize the memory of each Crow 
4: While (iters itermax) do 

ere for i= 1: to Npdo 
6 . 
7 
8 





Randomly choose one of the crows to follow (for example, j) 
Define an awareness probability 
if (r,>AP /"") then 


9: xbitertloy Liter, Yi afpitter a(mbtter_y biter ) 

10: else 

11: x*ier+l—a random position of search space. 
12: end if 

13: end for 

14: Check the feasibility of new positions 

152 Evaluate the new position of the Crows 
16: Update the memory of crows 

Af: end while 





4. RESULTS AND DISCUSSION 

In this section, we present the results obtained based on CSA for solving the ED problem and 
compare this results with the CM (Conventional Method) [15] and Particle Swarm Optimization [32]. A six 
units power unit system to explore our idea on using CSA to find the optimal set of power generation of the 
system. CSA will be used in this paper to solve the problem of economic dispatch of a test network of 26 
nodes and the morocco network. The programs are developed in MATLAB 7.9 environment. 

The adopted system is expected to produce demand power of 700 MW. The tuning parameters for 
CSA are given in Table 1, the Table 2 shows the cost coefficient of the six generators, under study, while the 
matrix is the loss coefficient matrix of the six units power system. From the results of Table 3, we notice that 
CSA give us the same production cost, and CM gives a slightly lower cost of $ 0.08 / h, by cons, delta given 
by CM is larger than that given by CSA, mean it is less accurate although we also see that it is the fastest. 
CSA gives us a good production cost and good accuracy. 


Table 1. parameters of CSA [34] 
Algorithm/parameters AP fl 
CSA 0.2 2 











In Figure 1, we show the convergence of the metaheuristic search process based on CSA in both the 
best and average cases. To see the difference between our new approach and another known method, we will 
compare the production cost found by CSA to that found by PSO [33]. The comparison is represented by the 
graph in figure 2. It can be seen that CSA provided the minimum fuel cost in this case compared to other 
reported methods in the literature. This shows that the CSA is more effective in finding the best load for the 
three generator system. 

In this case, we will test the operation of CSA. For this, we will use a simple network of 26 nodes 
with 6 production units. The total demand of the network is equal to 700 MW and loss coefficients are as 
follows: 


14 1.7 15 1.9 26 2.2 
1.7 60 13 16 1.5 2.0 
5/15 13 65 1.7 24 1.9 
1.9 16 1.7 7.2 3.00 2.5 
26 15 24 3.0 69 3.2 
2.2 2.0 1.9 2.5 3.2 8.5 


Bij = 107 


The simulation results are presented in Table 3. 
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Table 2. The parameters of the cost function and generators limits of the six-unit system 








Units a ($/MW2) b ($/MW) c ($) Pmin (MW) Pmax (MW) 
1 0.15240 38.53 756.79886 10 125 
2 0.10587 46.15916 451.32513 10 150 
3 0.02803 40.39655 1049.9977 35 225 
4 0.03546 38.30443 1242.5311 35 210 
5 0.02111 36.32782 1658.5696 130 325 
6 0.01799 38.27041 1356.6592 125 315 





Table 3. Results of the economic dispatching of six-unit system 








CM PSO CSA 

Pl (MW) 29.4552 28.28 29.45 
P2 (MW) 10 10 10 
P3 (MW) 118.814 119.02 118.702 
P4 (MW) 118.420 118.79 118.339 
P5 (MW) 230.559 230.78 230.441 
P6 (MW) 212.511 212.56 212.394 
Pl (MW) 19.3361 19.4319 19.3128 
Delta (MW) 0.000702 0.006688 0.000677 
Fuel cost ($/h) 36926 36912.16 36906 
t (s) 64.61 43.72 43.45 
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Figure 1. The convergence of CSA in the case of economic dispatch of six-unit system 
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Figure 2. Comparison graph between CSA and PSO 


5. CONCLUSION 

In this paper, we proposed a crow search algorithm to solve the problem of economic dispatch of the 
morocco electricity network. The practicality of the proposed metaheuristics CSA was tested for six power 
generators test case. The gained results were compared to existing results based on PSO and CM methods. It 
was shown that CSA are superior in obtaining a combination of power loads that fulfill the problem 
constraints and minimize the total fuel cost. CSA found to be efficient in finding the optimal power 
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generation loads. CSA was capable of handling the non-linearity of ED problem. The evolved power using 
CSA minimized both the cost of generated power, the total power loss in the transmission and maximizes the 
reliability of the power provided to the customers. The programs were developed using MATLAB and tested 
a network of 26 nodes. The results have shown that our CSA to give us a better performance with optimal 
results in all cases and respecting the constraints imposed. 


ACKNOWLEDGEMENTS 
The authors are very much thankful to the unanimous reviewers of the paper and editors of the 
journal for their constructive and helpful comments that improved the quality of the paper. 


REFERENCES 

[1] C.He-Rui, P. Xu, "Study on Smart Grid System Based on System Dynamics", TELKOMNIKA Indonesian Journal 
of Electrical Engineering Vol. 12, No. 12, pp. 7979-7986, December 2014. 

[2] Shahinzadeh H, Hasanalizadeh-Khosroshahi A. "Implementation of Smart Metering Systems: Challenges and 
Solutions". TELKOMNIKA Indonesian Journal of Electrical Engineering. 2014; 12(7). 

[3] Hermagasantos Z, Yusra S, Ali M. Implementation of Electricity Competition Framework with Economic Dispatch 
Direct Method. Telkomnika. 2012; 10(4): 667-674. 

[4] Wollenberg B, Wood A (1996) Power generation, operation and control. Wiley, New York 

[5] Dieu VN, Schegner P, Ongsakul W (2013) Pseudo-gradient based particle swarm optimization method for 
nonconvex economic dispatch. In: Power, control and optimization. Springer, New York, pp 1-27. 

[6] Muhtazaruddin, M.N.B., Bani, N.A., Aris, S.A.M., Kaidi, H.M., Fatah, A.Y.A., Jamia, J.J., Muhammad-Sukki, F. 
and Abu-Bakar, S.H. “Distribution power loss minimization via distributed generation, capacitor and network 
reconfiguration”. Indonesian Journal of Electrical Engineering and Computer Science. 2017; 5(3): 488-495 

[7] Arriffin, A.M., Othman, M.M., Kamaruzaman, A.A.M., Musirin, I., Yahya, A., & Latip, M.F.A. “Stochastic 
Approach of Voltage Optimization to Maximize Power Saving in a Building”. Indonesian Journal of Electrical 
Engineering and Computer Science. 2017; 8(1): 268-272. 

[8] Moradi-Dalvand M., B. Mohammadi-Ivatloo, A. Najafi, A. Rabiee. Continuous quick group search optimizer for 
solving non-convex economic dispatch problems. Electric Power Systems Research. 2012, Vol 93, pp. 93-105. 

[9] L. Kirchmayer, Economic Operation of Power Systems. New Delhi: Wiley Eastern Limited, first ed.,1979. 

[10] A. J. Wood and B. F. Wollenberg, Power Generation Operation and Control. New York: John Wiley & Sons, Inc., 
second ed., 2006. 

[11] G. F. Reid and L. Hasdorff, “Economic dispatch using quadratic programming,” IEEE Transactions on Power 
Apparatus and Systems, vol. PAS-92, no. 6, pp. 2015— 2023, Nov 1973. 

[12] K. Aoki and T. Satoh, “New algorithms for classic economic load dispatch,” IEEE Transactions on Power Apparatus 
and Systems, vol. PAS-103, no. 6, pp. 1423— 1431, June 1984 

[13] J. K. Delson and S. M. Shahidehpour, “Linear programming applications to power system economics, planning and 
operations,” IEEE Transactions on Power Systems, vol. 7, no. 3, pp. 1155-1163, Aug 1992 

[14] S. Subramanian and S. Ganesan, “A simple approach for emission constrained economic dispatch problems,” 
International Journal of Computer Applications, vol. 8, no. 11, pp. 39-45, October 2010, published By Foundation 
of Computer Science. 

[15] D. D. Obioma and A. M. Izuchukwu, “Comparative analysis of techniques for economic dispatch of generated 
power with modified lambda-iteration method,” in Proceedings of the 2013 IEEE International Conference on 
Emerging Sustainable Technologies for Power ICT in a Developing Society (NIGERCON), Nov 2013, pp. 231-237 

[16] S. K. Mishra and S. K. Mishra, “A comparative study of solution of economic load dispatch problem in power 
systems in the environmental perspective,” Procedia Computer Science, vol. 48, pp. 96 —100, 2015. 

[17] Z. L. Gaing, “Particle swarm optimization to solving the economic dispatch considering the generator constraints,” 
IEEE Transactions on Power Systems, vol. 18, pp. 1187-1195, 2003. 

[18] H. Faris, A. Sheta, and E. “Oznergiz, “Modelling hot rolling manufacturing process using soft computing 
techniques,” International Journal of Computer Integrated Manufacturing, vol. 26, no. 8, pp. 762-771, 2013 

[19] A. F. Sheta, H. Faris, and E. “Oznergiz, “Improving production quality of a hot-rolling industrial process via genetic 
programming model,” Int. J. Comput. Appl. Technol., vol. 49, no. 3/4, pp. 239-250, Jun. 2014. 

[20] H. Faris and A. F. Sheta, “A comparison between parametric and non-parametric soft computing approaches to 
model the temperature of a metal cutting tool,” International Journal of Computer Integrated Manufacturing, vol. 29, 
no. 1, pp. 64-75, 2016. 

[21] M. Braik, A. Sheta, and A. Ayesh, “Particle swarm optimization enhancement approach for improving image 
quality,” Int. J. Innov. Comput. Appl., vol. 1, no. 2, pp. 138-145, Jan. 2007. 

[22] B. Solaiman and A. Sheta, “Evolving a hybrid KMeans clustering algorithm for wireless sensor network using PSO 
and GA,” International Journal of Computer Science Issues, vol. 12, no. 1, pp. 23-32, 2015. 

[23] J. Sasikala and R. M, “Optimal _ based economic emission dispatch using simulated annealing,” International 
Journal of Computer Applications, vol. 1, no. 10, pp. 55-63, February 2010, published By Foundation of Computer 
Science. 


Indonesian J Elec Eng & Comp Sci, Vol. 13, No. 1, January 2019 : 347 — 353 


Indonesian J Elec Eng & Comp Sci ISSN: 2502-4752 im) 353 





[24] H. Hardiansyah, “Artificial bee colony algorithm for economic load dispatch problem,” IAES International Journal 
of Artificial Intelligence, vol. 2, no. 2, pp. 90— 98, 2013 

[25] W.M. Aly and A. Sheta, “Parameter estimation of nonlinear systems using levy flight cuckoo search,” in 
Proceedings of the Thirty-third SGAI International Conference on Innovative Techniques and Applications of 
Artificial Intelligence, Cambridge, England, UK, 2013, pp. 443-449 

[26] “Comparative analysis of optimal load dispatch through evolutionary algorithms,” Ain Shams Engineering Journal, 
vol. 6, no. 1, pp. 107-120, 2015 

[27] Ajenikoko, G.A., & Olabode, O.E. “Optimal Power Flow with Reactive Power Compensation for Cost and Loss 
Minimization on Nigerian Power Grid System”. Indonesian Journal of Electrical Engineering anda Informatics 
(JEED). 2017; 5(3): 236-247. 

[28] Jiang A, Ertem S, Subir S, Kothari D (1995) Economic dispatch with non-monotonically increasing incremental cost 
units and transmission system losses. Discussion. IEEE Trans Power Syst 10:891—897 

[29] A. Askarzadeh, “A novel metaheuristic method for solving constrained engineering optimization problems: crow 
search algorithm,” Computers & Structures, vol. 169, pp. 1-12, 2016 

[30] A. Askarzadeh, “Electrical power generation by an optimised au-tonomous pv/wind/tidal/battery system,” IET 
Renewable Power Gener-ation, 2016. 

[31] D. Oliva, S. Hinojosa, E. Cuevas, G. Pajares, O. Avalos, and J. Galvez,’ “Cross entropy based thresholding for 
magnetic resonance brain images using crow search algorithm,” Expert Systems with Applications, vol. 79, pp. 164— 
180, 2017 

[32] A. L. Devi and O. V. Krishna, “Combined economic and emission dispatch using evolutionary algorithms: A case 
study,” ARPN Journal of Engineering and Applied Sciences, vol. 3, no. 6, pp. 28-35, 2008 

[33] GAING, Z. L. Particle Swarm Optimization to Solving the Economic Dispatch Considering the Generator 
Constraints. IEEE transactions on power systems, 2003, vol. 18 (3), pp. 1187-1195 

[34] R . Habachi et al .Eagle Strategy Based Crow Search Algorithm for Solving Unit Commitment Problem in Smart 
Grid System Indonesian Journal of Electrical Engineering and Computer Science Vol. 12, No. 1, , pp. 17~29 
October 2018 





Resolution of economic dispatch problem of the moroccan network ... (Rachid Habachi) 


